Corrections 


Page  7,  line  3,  should  read  as  follows: 
in  terms  of  elements  of 

Page  13,  line  4,  should  read  as  follows: 


Page  16,  line  8,  should  read  as  follows: 


Page  19,  line  1,  should  read  as  follows: 

^  Aft,*  "  '  W\-io 

Page  22,  line  7,  should  read  as  follows: 

at  depth  £  at  the  bottom  of  the  layer  is  4k 


Page  22,  line  10,  should  read  as  follows: 

Page  28,  line  2,  should  read  as  follows:  but  from  the  relation 

A A  - *X  f^w  A3|  -  i 
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ABSTRACT 

This  report  develops  by  matrix  methods  the  integral  solution 
of  the  wave  equation  for  point  sources  of  harmonic  waves  in  a 
liquid  layer  of  a  multilayered  liquid-solid  half  space  in  a  form 
convenient  for  numerical  computation  on  a  high-speed  digital 
computer.  Only  the  case  is  considered  of  a  high-speed  liquid 
bottom  underlying  the  stack  of  layers.  The  integral  over  wave 
number  has  singularities  in  the  integrand  and  is  conveniently 
transformed  into  the  complex  plane.  By  a  proper  choice  of 
contours,  complex  poles  are  displaced  to  an  unused  sheet  of  the 
two-leaved  Riemann  surface,  and  the  integral  solution  for  the 
multilayered  system  reduces  to  a  sum  of  normal  modes  plus  the 
sum  of  two  integrals,  one  along  the  real  axis  and  the  other 
along  the  imaginary  axis.  Both  integrals  are  evaluated  by  a 
Gaussian,  quadrature  formula.  Sample  computations  are  presented 
for  low-frequency  propagation  in  the  Arctic  Ocean  sound  channel. 
These  are  preliminary  computations  and  the  ice  layer,  which 
averages  three  meters  in  thickness,  is  not  included  in  the 
layered  system.  The  effects  of  the  ice  layer  on  propagation 
are  currently  under  investigation. 


FIGURES 


Figure  1. 
Figure  2. 
Figure  3. 

Figure  4. 

Figure  5. 


Figure  6. 


Figure  7. 

Figure  8. 
Figure  9. 
Figure  10. 


Multilayered  half-space. 

Contours  for  integration. 

Variation  of  sound  velocity  with  dej  .h  for  Model  A. 

Table  2  gives  additional  parameters  for  this  model. 
Computations  for  Model  A.  Phase  -  and  group-velocity 
dispersion  and  excitation  function  of  pressure 
dependent  only  on  layering. 

Computations  for  Model  A.  Variation  with  range  of  the 
absolute  value  of  pressure  for  the  normal-mode  contri¬ 
bution  of  the  sound  field.  Source  depth  150  m.  Hydro¬ 
phone  depth  50  m.  Source  frequency  10  Hz.  Source 
pressure  amplitude  1  ayne/cm^  re  Ini. 

Computations  for  Model  A.  Variation  with  range  of 
absolute  value  of  pressure  for  the  integral  contribution 
of  the  sound  field.  Source  and  detector  same  as  for 
Figure  5. 

Computations  for  Model  A.  Variation  with  range  of  absolute 
value  of  pressure  of  the  total  sound  field.  Source  and 
detector  same  as  for  Figure  5. 

Same  as  Figure  5  but  computations  carried  to  longer  range. 

Same  as  Figure  6  but  computations  carried  to  longer  range. 

Same  as  Figure  7  but  computations  carried  to  longer  range. 
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DEFINITION  OF  SYMBOLS 

compressional-wave  velocity  in  the  m-th  layer 

shear-wave  velocity  in  the  m-th  layer 

thickness  of  the  m-th  layer 

vertical  coordinate 

range  between  source  and  detector 

time 

angular  frequency 
phase  velocity 
wave  number 

Nh1 

Imaginary  part  of  a  complex  number  x  +  iy 
velocity  potential  in  the  m-th  layer 
normal  stress  in  the  m-th  layer  parallel  to  z  axis 
pressure  in  the  m-th  liquid  layer 
horizontal  particle  displacement  in  the  m-th  laye*' 
vertical  particle  displacement  in  the  m-th  layer 
horizontal  particle  velocity  in  the  m-th  layer 
vertical  particle  velocity  in  the  m-th  layer 
density  in  the  m-th  layer 
density  at  the  source. 


Jo 

X 

K0 

u. 


Eessel  function  of  order  0 
Y  bessel  function  of  order  0 
K  Bessel  function  of  order  0 
Hankel  function  of  the  first  kind  of  order  0 
SO;  Hankel  function  of  the  second  kind  of  order  0 


(Yv\ 


fV*  ^ 


kA rv\^ 


“°^AYV. 


St 


-  k*  s  V  <  V* 


>k>  v 


/wv. 


c  =  • 

Arvv\  ^ 


1 


\  -/  Sl. 


/TT1 

vP^ryvv. )  ^ 


C  <  4 


/yA» 
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^  -  -  <■' 


\  -IS. 
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JsySB*«oW{^»A  ✓  "** 
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5. 


INTRODUCTION 


This  report  develops  the  integral  solution  of  the  wave 
equation  for  point  sources  of  harmonic  waves  in  a  liquid 
layer  of  a  multilayered  liquid-solid  half  space  in  a  form 
convenient  for  numerical  computation  on  a  high-speed  digital 
computer.  Only  the  case  of  a  liquid  half- space  underlying 
the  stack  of  layers  is  considered.  Furthermore,  it  is  assumed 
that  the  speed  of  sound  in  this  last  layer  of  infinite  thick¬ 
ness  is  greater  than  the  speed  of  sound  in  liquid  layers 
immediately  above  it.  The  problem  of  Immediate  concern  is 
low-frequency  propagation  from  harmonic  sources  in  the  central 
Arctic  Ocean  at  ranges  up  to  ten  times  the  water  depth  in 
shallow  and  Intermediate  water  depths.  Solutions  are  develop¬ 
ed  for  the  pressure  perturbations  detected  by  a  hydrophone  at 
depth  or  the  ice  vibrations  detected  by  a  seismometer  on  the 
ice  surface.  A  modification  of  the  formulation  Is  possible  to 
include  harmonic  vibration  sources  on  or  in  the  Ice,  although 
it  appears  at  the  present  time  that  surface  sources  are  of 
limited  application  (Greene,  1968). 

The  solution  of  the  wave  equation  presented  here,  based 
on  the  Thomson-Haskell  matrix  method  (Thomson,  195G;  Haskell, 

1953),  follows  Harkrider  (1964)  for  the  solution  of  the  wave 
equation  in  an  n-layered  solid  half  space.  Layer  matrices  of 
the  type  given  by  Dorman  (1962)  for  computing  dispersion  in 
an  n-layered  liquid  -  solid  half-space  are  used  for  the 
liquid  layers.  An  application  of  the  theorem  that  the  inverse  of  the 
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oroduct  matrix  for  the  layered  system  above  the  source  has  the  same 
form  as  the  inverse  of  a  layer  matrix  reduces  the  integrand 
of  the  integral  solution  to  a  simple  form  in  terms  of  product 
matrices  derived  in  the  source-free,  plane-wave  case.  The 
integral  over  wave  number  from  zero  to  infinity  has  singular¬ 
ities  in  the  integrand  and  is,  therefore,  conveniently  trans¬ 
formed  into  the  complex  t,  -  k  +  It  plane.  Performing  the 
contour  Integration  in  the  ;  plane  and  noting  that  branch 
line  integrals  corresponding  to  branch  points  generated  by 
each  layer  matrix  except  the  last  are  zero,  we  are  left  with 
an  expression  analogous  to  that  of  Sorensen  (1955)  and  Leslie 
and  Sorensen  (1961)  for  the  two-layer  liquid  half-space  with 
a  high-speed  bottom.  By  the  proper  choice  of  contours, 
complex  poles  are  displaced  to  an  unused  sheet  of  the  two-leaved 
Riemann  surface,  and  the  integral  solution  for  the  multilayered 
system  reduces  to  a  summation  of  normal  modes  plus  the  sum  of 
two  integrals  over  wave  number,  one  along  the  real  z,  axis  and 
the  other  along  the  imaginary  5  axis.  Both  of  these  integrals 
are  conveniently  computed  by  a  Gaussian  quadrature  formula. 

The  physical  interpretation  of  the  final  solution  is 
straightforward.  The  normal-mode  terms  correspond  to  waves 
trapped  in  the  Arctic  sound  channel;  that  is,  refracted  and 
surface  reflected  (RSR)  sounds  and  reflected  sounds  incident 
on  the  bottom  beyond  the  critical  angle.  The  Integral  over 
wave  number  along  the  real  axis  corresponds  to  waves  incident 


Q 

-*■  • 

on  the  bottom  at  angles  greater  than  the  critical  angle  of  total 
reflection.  Since  energy  leaves  the  guide  continuously  as  waves 
travel  down  the  guide,  this  term  is  of  importance  only  to  ranges 
of  about  ten  times  the  water  depth.  The  contribution  to  the  sound 
field  of  the  integral  along  the  imaginary  axis  is  of  importance  only 
very  near  the  source  since  the  integrand  decays  exponentially  with 
range  and  wave  number. 

The  usefulness  of  the  normal  mode  terms  for  describing  long- 
range  explosive  sound  transmission  in  the  central  Arctic  Ocean 
is  shown  by  Kutschale  (1969) .  In  that  work  the  present  formulas 
for  the  normal  modes  were  extended  to  explosive  sources  and  the 
Fourier  integral  for  each  mode  was  evaluated  by  the  principle  of 
stationary  phase  (Peker.is,  1948)  .  Attenuation  by  the  rough  ice 
boundaries  was  also  included,  although  omitted  here. 

FORMAL  SOLUTION 


Source  Free  Case 

Consider  the  r-layered  interbedded  liquid-solid  half-space 
shown  in  Figure  1  in  which  the  last  layer  of  infinite  thickness  is 
liquid  and  has  a  higher  sound  velocity  than  liquid  layers 
immediately  above  it;  that  is,  a  high-speed  bottom.  A  point 
source  of  harmonic  waves  is  located  in  one  of  the  liquid  layers. 
The  velocity  potential  in  the  m-th  layer,  which  is 
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where  A  and  A^  are  constants. 


Placing  the  origin  at  the  (m-l)-st  Interface  we  get 


/\  I 

M  =  LUjf'1^  C  AivvJ') 


for  z  =  0 


^ftA.  ~-i  k*'  (  Afy^-V  Aftrv^ 

(JhA^"  ( V* k*) cos^v 

—  ^Pav\  (,  A|V\“  Afwv')  Ss'  n.  1  Arv 

for  ^ 


^  . 


Substituting  expressions  for  [J^  -  A^]  and  [A^  +  A^]  from 
equations  (1)  in  equations  (2)  yields 


13. 


since  the  vertical  particle  velocity  and  pressure  are  continuous 


across  each  interface.  For  the  n-th  layer 


A.„+a1  -  Tii’Vi  _  A avur( 


lk, 


or  in  matrix  notation 


C  cO  pv-\ 


uj 
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An  + A* 
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o  -L 

Applying  boundary  conditions  that  the  pressure  vanishes 
at  the  surface  and  that  no  upgoing  waves  travel  from 
Infinity,  (pzz)0  =  0,  =  0,  the  relation  (3)  reduces  to 


»  A.  —v 

A  ' 

«■ 

-  E  A 

o 

Consider  now  a  solid  layer  between  two  liquid  layers  or  at 
the  surface  of  the  laminated  halfspace.  For  this  layer  we  may 

write  (Thomson,  1950;  Dorman,  1952)  , 
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15. 


tvhere  the  fourth  element  of  the  column  vector  is  the  tangential 
stress  which  is  zero  at  a  solid-liquid  boundary.  The  matrix 
elements  are  given  in  the  appendix  and  may  be  derived  following 
Haskell  (1953).  Equation  (4)  yields  the  three  equations. 


u. 


-  ur.  _ 
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<WC 


*  A.  * 


(6) 

^  » 1 


Substituting^^  from  equation  (5)  in  equations  (6)  and  (7) 
yields 


which  is  the  same  as  the  matrix  relation  for  a  liquid  layer. 
In  general,  then,  for  interbedded  solids  and  liquids 
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Point  Source  of  Harmonic  Waves  In  a  Liquid  Layer . 

Divide  the  source  layer  into  two  layers  as  shown  In 
Figure  1.  At  £*D  the  pressure  is  continuous.  The  vertical 
particle  velocity  is  continuous  everywhere  in  the  plane  defined 
by  z  =D  except  at  the  point  source  where  the  liquid  above  and 
below  the  source  moves  in  opposite  directions.  This  may  be 
expressed  by  writing 

-  3.  K  . 


In  matrix  notation  for  the  source  layer 


For  the  layers  below  z  =D  we  have 


and  for  the  layers  above  z=D 


18. 


where  A<~a=  '  Q-g^  and  A.-  CL^  Q-j  •  •  •Q_v 


Returning  now  to  the  relation  for  the  n-th  liquid  layer 
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X-  £r0 


ft 


We  now  define 


Hence  (12)  may  be  written  as 

L*0  nr  k,x*y] 

-A^-  cr„x  *  Ju't 

a;  -  Ta,X  +  T,^  . 


Adding  equations  ( 1 A )  we  get 

0  =  (J1]L  +  J21)  X  +  (J12  +  J22)  Y  and  if  we 
=  J-q  +  J  and  V  =  J  +  and  solve  fcr  X  we  hav? 


let  T 
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22. 


Therefore  the  integral  solution  for  the  surface  verticle  particle 
velocity  is  UJ"0  (  ^  ^  t)  - 
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so 
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and  from  the  formula 

A  % 
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Evaluation  of  the  Integral  Solutions. 


The  Integral  solutions  have  singularities  corresponding 
to  zeros  of  A^\ 


4- 


i*  -  C  *)pn 

These  are  simple  poles  for  complex  wave  number  and ,  therefore, 
the  Integral  solutions  are  conveniently  evaluated  by  contour 
integration  in  the  complex  ^  =  k  +  i ^  plane.  The  contours 

are  shown  in  Figure  2.  First  the  transformation 


lA\ 


In  addition  to  simple  poles  at 


is  made. 


An 


A** 


branch  points  occur  at 


vlr  * 


-  o 


for  the  liquid  layers  and 


=  o  -  O 

for  each  solid  layer.  But  since  the  integrands  are  even 

J 


functions  of 


i  <r  y 


7 


and 


branch  line  integrals  corresponding  to  branch  cuts  made  to 
these  branch  points  cancel  and  only  the  integral  corresponding 
to  .the  branch  cut  for  -  O  must  be 

considered.  »  ^ 


p  J  &ra**»4*r>a**«*-;‘* 


25. 


The  Reimann  surface  h^s  tv/o  sheets.  To  satisfy  the  varishinn 


of  ULT0  as 


^  approaches  infinity,  we  nust  remain  on  the  sheet  of  the 


Rienann  surface  v/herr  the  real  part  of  L 


is  greater  than  zero. 


For  tie  contours  shovm  in  Figure  2  complex  poles  of 

Aw _  +  __ 


have  been  displaced  to  an  unused  Relmann  sheet  and  all  poles 

are  on  the  real  k-axis.  See  Hwino,  Jardetzkv  and  Fress  (1957  , 
pages  135  to  137)  for  the  proof  for  a  two -layer 

liquid  half  space.  Along  the  branch  cut 


o 


Is  pure  imaginary  and  along  the  real  axis  from 


infinity 


In  the  first  quadrant 


is  real. 

is  positive  and 


in  the  fourth  quadrant  negative.  The  integral  for 

,  ,®  CCdt 

(omitting  the  term)  is  now  transformed  to 


where 


in  the  first  quadrant,  (17)  is 


/°° 


Likewise,  since 


y° 

(  X  t  =  o 

\to  /° 

•  o  .  tfl0 

approaches  zero  as  r  approaches 


infinity  in  the  fourth  quadrant  (18)  is 

y-  LOO  o 


o  co  o 


<3TT(lJ'Res(Xa) 
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In  the  first  integral 
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and  therefore 
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but  ^2.1  ^22  "*  A^2  ^21  c  ^ 
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Likewise,  in  the  second  integral 


Ft<f^>FKJk5i?V 
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Therefore,  the  integral  solution  reduces  to 


29. 


From  the  formulas 


^Ko(.t0  and 


»  T.lVei  -  c  V-f 

and  multiplying  jjJ"  by  '  O  for  a  constant  pressure 

0  LLOf^ 


source  of  pressure  amplif.de  PQ  at  unit  distance  from  the  source 

W7’0""  \ 


Q.  (><0  Udk. 


In  summarv,  the  expressions  for  vertical  particle  velocity 
at  the  surface,  vertical  particle  displacement  at  the  surface, 
and  pressure  in  a  liquid  layer  are: 

#“ 

f  f - \?(A  ^ \T0 -  C>(  jkcS^kcfr 


has  been  used. 


33. 


For  programming  It  Is  convenient  to  write: 


wr0— ' W,+Va-cV/j+cV/4 


and  for  the  absolute 


value  of  vertical  particle  velocity,  a  quantity  conveniently 
measui*ed  in  transmission  experiments, 

Kol-  ( C-H  *v/*T  *  ( 
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^  a?0  p.f^X^\H^t\^TD(kQkdk 
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3JJ. 


'  2c\X^(c^  pH  . 


w=M| 

k:«C( 


Likewise  for  vertical  particle  displacement 


The  layer  matrices  for  the  integral  along  the  real  axis 
are  programmed  as 


-Xfw^<XjQ 

^\\  -X4a^  IP^3  X^^] 

HKl  Xfw^*^}  (0*\( 

for  solid  layers  and 

[a^ 

-  w 

for  liquid  layers.  For  the  integral  along  the  imaginary  axis 


the  matrices  are 


for  liquid  layers.  Likewise,  for  the  ndrmal  mode  terms  the 
layer  matrices  are  programmed  as 


39. 
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X^q^ 
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T^Ari 
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L  IH, 

4h^ 
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for  solid  layers  and 

tPOAa 

j 

VQ-<w^u 

for  liquid  layers. 


40. 


NUMERICAL  COMPUTATIONS 


Computer  programs  were  written  in  double  precision  Fortran 

.  2 

IV  to  compute  the  pressure  in  dynes/cm  and  the  vertical  particle 
displacement  in  millimicrons.  The  programs  are  run  on  the  IBM 
360/91.  For  numerical  results  presented  here,  computing  time 
was  under  six  minutes.  It  appears  that  practical  limits  of 
numerical  precision  make  the  present  develop  ':t  most  useful  for 
the  Arctic  at  freeruencies  below  50  Hz  in  water  depths  to  1  km. 

At  higher  frequencies  or  in  greater  water  depths  the  integrands 
may  be  so  oscillatory  that  it  is  difficult  to  achieve  the  desired 
accuracy  in  the  numerical  integrations  without  excessive  computing 
time. 


Computations  are  made  in  two  stages.  The  first  program,  an 
extension  of  Dorman’s  (1962)  PV  7  dispersion  program,  computes 
phase-  and  group-velocity  dispersion,  the  excitation  function  depen¬ 
dent  only  on  layering  of  the  medium,  and  the  excitation  function  for 
the  particular  source  and  detector  depths.  In  the  first  case  for 
the  m-th  normal  mode  the  excitation  function  is  defined  by 


Wv 
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for  pressure  and 

C  36 ^{uf 
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for  vertical  particle  velocity.  In  the  second  case  the  excitation 


function  is  —  *— «.  ►k  —  _ _  ^ 


ps  C/a  A„  1  \  C^A"  _  +  ‘W  Aa^\l 

1  ? «  J 


for  pressure  and 


9-  ju?  T  ^  (A  ^  ^ 


p^p&C^  A„  ^  I  C"A" 

3c  L;Twir 


u^  .  0«H 

C”  Cmi 


^  ^UX<r0  ^  ^w, 

<^c* 

for  vertical  particle  velocity.  These  definitions  were  chosen  tc  be 
useful  also  at  long  ranges  where  uriko  = 

w  ^ 

The  second  program  computes  the 
three  integrals,  the  normal  mode  contribution  to  the  sound  field, 


and  writes  and  plots  the  absolute  value  of 


i  ' 


or  wQ  as  a 


I 
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Computations:  Source 

*> 

**:  Source  Depth,  15.2  m; 
m 

Shear  Densi 

Velocity,  m/sec  gm/cm 

0  1.03 


0 


1.24 


'pi'll  WPI^HPiUPIJ  l  JU !!'i!P 


TABLE  2 
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Integral  Real  Axis,  Real  Part 


Frequency, 

Hz 

Sorensen's  (1959)  Value, 

dynes/cm^ 

Our  Value, 
dynes/cm2 

10 

,0014713 

.004685 

20 

.030170 

.030153 

40 

.040278 

.040485 

80 

.023193 

.023121 

160 

-  .021839 

-  .022201 

320 

-  .074395 

-  .074011 

;  » 

-.1 


I  -5 
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TABLE  3. 

Integral  Real  Axis,  Imaginary  Part 


Frequency, 

Hz 

Sorensen's  Value, 
dynes/cm2 

Our  Value, 
dynes/cm2 

10 

-  .002567 

-  .002551 

20 

-  .022260 

-  .0222*10 

40 

-  .061960 

-  .062274 

80 

-  .050986 

-  .050641 

160 

-  .0*17979 

-  .047750 

320 

.006299 

.008316 

TABLE  4. 


Integral  Imaginary  Axis 


Frequency, 

Hz 

Sorensen’s  Value, 
dynes/cm^ 

Our  Value 
dynes/cm2 

10 

.095080 

.095146 

20 

.081938 

.082076 

40 

.041887 

.041759 

80 

.030795 

.030937 

160 

.030214 

.030356 

320 

.010593 

.011043 

TABLE  5. 
Model  A 


Layer 

Thickness,  m 

Compressional 
Velocity,  m/sec 

Shear  Velocity, 
m/sec 

Density, 

gic/cm3 

50 

lb3k. 0 

0 

1.03 

50 

1K37.0 

0 

1.03 

50 

lMO.O 

0 

1.03 

50 

li»i»3.o 

0 

1.03 

50 

lUl|6.0 

0 

1.03 

50 

1U50.0 

0 

1.03 

50 

1^53.2 

0 

1.03 

Infinite 

1600.0 

0 

1.20 

FREQUENCY 


HANOE  (M) 
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function  of  range.  In  this  program  provision  is  also  made  tc 
compute  the  three  integrals  as  a  function  of  frequency  or  detector 
depth  at  a  specified  range.  The  integrals  are  evaluated  by 
Gaussian  quadrature  formulas  (see,  for  example,  Hildebrand,  1956). 
Care  must  be  taken  to  employ  a  sufficient  number  of  Gaussian  points 
to  obtain  the  desired  accuracy  of  the  integrals.  A  thirty- two 
point  Gaussian  quadrature  formula  was  found  to  give  sufficient 
accuracy  for  the  computations  presented  here  out  to  a  range  of 
1500  m,  but  at  longer  ranges,  higher  frequencies,  or  greater  water 
depths  a  ninety-six  point  formula  is  used.  Abscissas  and  weight 
factors  for  the  Gaussian  integration  are  given  by  Davis  and 
Rabinowitz  (1956,  1958).  The  integration  along  the  imaginary 
axis  extends  to  infinity,  but  in  practice  a  finite  upper  limit  is 
chosen  which  gives  sufficient  accuracy.  This  may  be  done  because 
the  integrand  is  an  exponential  function  of  -kr  for  large  k  or  r. 

The  programs  were  checked  by  computations  for  the  same  t\*o-layer 
model  of  Table  1  used  by  Sorensen  (1959).  Tables  2,  3  and  4  show 
that  the  two  sets  of  numerical  results  are  in  very  close  agreement. 

As  an  illustration  of  the  progrms  for  the  Arctic,  computations 
were  made  for  the  layered  model  A  of  Figure  which  closely  follows 
the  observed  variation  of  sound  velocity  with  depth.  Additional 
parameters  for  the  model  are  given  in  Table  5.  This  is  a  prelim¬ 
inary  model  and  the  ice  sheet ,  which  averages  about  3-m  in  thick¬ 
ness  in  the  central  Arctic  Ocean,  was  not  included.  The  effects  of 
the  ice  sheet  on  propagation  are  currently  under  investigation. 
Figure  4  shows  phase  -  and  group-velocity  dispersion  and  the 
excitation  function  for  pressure.  At  a  frequency  of  10  Hz  only 
two  normal  modes  are  excited.  The  range  dependence  of  the  absolute 
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value  of  pressure  for  the  normal  mode  contribution,  the  Integral 
contribution,  and  the  total  pressure  are  shown  In  Figures  5  through 
10. 
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APPENDIX 

Matrix  elements  for  normal  modes  (see  Haskell,  1953,  and 
Dorman,  1962). 

Solid  layers 

(T  W'),^  ~  (A.  ivn)^  =  “ (p/v»\C-  ^  C  ^  <*?*  —  COSQitov') 

tA*»V  (fOn* * 

(,V~ 

^ArtwVa  + 

VC<*>Q* 

Ki,3  *  £ 1  f^c-1  H + 9^*^) 


P'>yA  Vi.'f/wr'XCosV/wt  -  Co^q„  ^ 

^A\»  = 

v- Ar  <“P^c’a-K  ^^wNT^-vCr^vif )%wfc\ 

Liquid  layers  «J 

^<w\\v<Q«»v\^=  C  os?** 

~  ”l^A  S \  Y\ 

__  p/w\C^* 

*  Vj 

Of  Am 

Matrix  elements  for  integrals  along  real  axis: 


Solid  layers 


k.^'amr, 


/^ik. 

V  V5^/wx  )  ' 


\<  *a \v\ 


rw 


^^a-i  '^^0,^=  ^-i^CoSV5,, 


+  \i5co^ ; 


{WU> 


»--L,a( 


i^v^V--s-ef(^C\Lv 

LtosV  -  “5>Q„1 


^lk  K^mUn 
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Liquid  layers 


<^S,x8, 

Vf 

(  afe\*~^fc*OTQ»« 

V  \Jsu*»)  J  rrr  -y 

C-OiU,  'k*«‘v 

t^V-  iSlm?, 

.  40 

va,*\  ^  iB^  %iv^l 

Matrix  elements  for  integral  along  imaginary  axis: 

Solid  layers 

CO^Q^ 


•^t. 


K*^ 


<vj*?  S«Q 


/W\ 


H'H,=  jcoVPM  -  C.OSQJ1 
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u  '  */W\ 


Liquid  layers 

~  COS"?/VW\ 

l£v\x=  <5>m^ 

°°fs 


C  P/vv\  Ul  S  \  \A 


Security  Classification 


DOCUMENT  CONTROL  DATA  •  R  &  D 


fS  re  uni/  s»f»raf»on  at  title,  body  cl  abstract  and  indexing  annotation  moif  fee  entered  when  the  overall  report  ia  claisilled) 


1  CiNAT’NY.  ACTIVITY  fCorporstt  author ) 

Lament -Doherty  Geoloyical  Observatory 
of  Columbia  University 

2ji.  REPORT  SECURITY  C  L  A  SSI  F  1C  A  7  iC'« 

Unclassified 

26.  GROUP 

Arctic 

"■Wfe '  integral  Solution  of  the  Sound  Field  in  a  Multilayered  Liquid-Solid 
Half-Space  with  Numerical  Computations  for  Low-Frequency  Propagation 
in  the  Arctic  Ocean. 

4  DElCRiPTivE  NOTES  (Type  ol  report  end  IncluMve  detea) 

Technical  Report 

5  *uTtiOP;5ifflMlni/nf,  middle  initial.  leet  nmtne) 

Henry  W.  Kutschale 

6  KURORT  DATE 

February  1970 

7a.  TOTAL  NO.  OF  PAGES  76.  NO.  OF  REFS 

64  13 

Be.  CONTRACT  OR  GRANT  NO 

N00014-67-A-0108-0916 

h  P  BO  i?r  t  N  C  . 

NR  307-320/1-6-69  (415) 

c. 

d. 

9«.  ORiGiNATOR'S  REPORT  NUMBER(S) 

1 

96.  other  REPORT  nOIS)  (Any  other  number*  diet  mmv  be  eaalgrted 
thie  report) 

DISTRIBUTION  STATEMENT 

Reproduction  of  this  document  in  whole  or  in  part  is  permitted  for 
any  purpose  of  the  U.S.  Government. 

11  SUPPLEMENTARY  NOTES 

12  SPONSORING  MILITARY  ACTIVITY 

U.S.  Naval  Ordnance  Laboratory, 

White  Oak,  Silver  Spring,  Maryland 
Office  of  Naval  Research. Washinaton. 

13  ABSTRACT 


This  report  develops  by  matrix  methods  the  integral  solution 
of  the  wave  equation  for  point  sources  of  harmonic  waves  in  a  liquid 
layer  of  a  multilayered  liquid-solid  half  space  in  a  form  convenient 
for  numerical  computation  on  a  high-speed  digital  computer.  Only  the 
case  is  considered  of  a  high-speed  liquid  bottom  underlying  the  stack 
of  layers.  The  integral  over  wave  number  has  singularities  in  the 
integrand  and  is  conveniently  transformed  into  the  complex  plane.  By 
a  proper  choice  of  contours,  comp] ex  poles  are  displaced  to  an  unused 
sheet  of  the  two-leaved  Riemann  surface,  and  the  integral  solution 
for  the  multilayered  system  reduces  to  a  sum  of  normal  modes  plus  the 
sum  of  two  integrals,  one  along  the  real  axis  and  the  other  along  the 
imaginary  axis.  Both  integrals  are  evaluated  by  a  Gaussian  quadrature 
formula.  Sample  computations  are  presented  for  low-frequency  propaga¬ 
tion  in  the  Arctic  Ocean  sound  channel.  These  are  preliminary  compu¬ 
tations  and  the  ice  layer,  which  averages  three  meters  in  thickness, 
is  not  included  in  the  layered  system.  The  effects  of  the  ice  layer  on 
propagation  are  currently  under  investigation. 
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